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摘 要 : 人 台风 相似 性 度量 方法 的 研究 对 防 灾 减 灾 、 辅 助 决策 等 具有 重要 意义 。 目 前 ， 人 台风 相似 性 的 研究 大 多 集中 在 台 
风 路 径 的 相似 性 度量 上 。 首 先 , 梳理 影响 台风 相似 性 度量 的 多 个 要 素 , 提出 了 基于 多 元 时 间 序 列 的 台风 数据 描述 方法 ; 


其 次 ， 提 出 了 台风 时 间 序 列 完整 性 、 一 致 性 评估 与 修复 方法 ; 最 后 ， 针 对 台风 时 间 序 列 的 不 等 长 问题 ， 
于 主 成 分 分 析 和 动态 时 间 弯 曲 距离 的 相似 性 度量 方法 。 通 过 实验 验证 ， 
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Abstract: The research on similarity measure of typhoon is significant to disaster prevention and auxiliary decision-making. At 
present, the similarity researches on typhoon mostly focus on the path. First, this paper combed the elements of similarity 
measure of typhoon and using multivariate time series to describe the typhoon due to the space-time and multi-factor 
characteristics of it. Second, it gave one method to measure and repair integrity and consistency of typhoon time series. Finally, 
in the view of unequal length of typhoon time series, this paper designed a method of similarity measure based on Principal 


Component Analysis and weighted dynamic time warping distance. Through the experiment, this method can realize effective 


0 ”引言 据 是 一 类 特殊 的 关系 数据 ， 具 有 时 间 属 性 ， 因 此 这 些 方法 不 适 
用 于 台风 数据 。 海 洋 领域 ， 文 献 [7-9] 通 过 抽样 方法 来 检验 海洋 
台风 作为 影响 我 国 最 大 的 海洋 灾害 之 一 , 仅 2016 年 , 共 造 数据 质量 是 否 合格 ， 并 没有 提出 修复 方法 ， 目 前 ， 对 海洋 数据 
成 174 人 死亡 、24 人 失踪 , 直接 经 济 损失 766.5 亿 元 叫 , 因此 ， 质量 评估 与 修复 研究 还 很 少 ， 而 海洋 数据 的 质量 高 低 直 接 影响 
研究 台风 对 防 灾 减 灾 、 辅 助 决策 等 具有 重要 意义 。 台 风 相 似 性 海洋 预报 等 相关 研究 。 因 而 ， 数 据 的 广泛 应 用 对 数据 质量 的 保 
度量 是 进行 路 径 预 报 、 灾 害 预 报 等 重要 手段 之 一 ， 研 究 其 可 以 障 提出 了 迫切 需求 。 
减少 人 员 伤亡 和 经 济 损 失 。 但 是 台风 数据 中 普遍 存在 质量 问题 ， 国内 外 研究 工作 者 对 时 间 序列 进行 深入 挖掘 ， 针 对 不 同 太 
导致 分 析 结 果 存 在 偏差 ， 无 法 提供 准确 的 信息 ， 造 成 灾害 预报 究 重点 提出 各 种 行 之 有 效 的 相似 性 度量 方法 , 大 致 可 分 为 两 类 ， 
普 误 ， 带 来 不 必要 的 财产 损失 和 人 员 伤 亡 ， 因 此 ， 评 估 人 台风 数 一 类 是 适合 等 长 时 间 序 列 ， 主 要 包括 模式 距离 、 余 弦 距离 和 欧 
居 质 量 并 对 其 进行 修复 是 相似 + 的 重要 一 部 分 。 式 距离 。 模 式 距离 0 是 表示 两 个 序列 趋势 的 差异 程度 ， 物 理 概 
台风 数据 质量 主要 体现 在 完整 性 和 一 致 性 两 个 方面 。 文 献 念 明 确 ， 划 分 合理 ， 但 表示 方法 粗糙 ， 结 论 不 够 精确 ; 基于 夹 
[2] 综 述 了 数据 完整 性 和 一 致 性 的 研究 现状 , 这 些 方法 主要 围绕 角 余 弦 距 离 中 计算 简单 ,但 匹配 序列 长 度 必须 相同 ; Agrawal 等 
医学 领域 ， 几 乎 没有 涉及 海洋 领域 。 李 建 中 [3-6] 团 队 在 数据 完 人 [9 提出 的 欧式 距离 (ED ), 计算 简单 , 时 间 复 杂 度 是 线性 的 ， 
整 性 和 一 致 性 方面 的 研究 大 多 针对 普通 的 关系 数据 ， 而 台风 数 应 用 较 多 ， 但 只 适用 于 等 长 的 时 间 序列 的 相似 性 度量 。 另 一 类 
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是 适合 
DTWU3] 是 
匹配 映射 的 相 化 
法 , 后 来 Berndt 等 人 0 将 它 用 于 时 间 序 列 的 相似 性 度 
的 基础 上 提出 了 最 大 特 
口 法 , 进行 心 电 信号 的 相似 性 度 
里 的 水 文 时 间 序 列 相似 性 
的 效率 ; 文献 [17] 提 出 
似 性 度量 方法 ， 较 好 地 解决 股市 技术 分 析 中 
出 了 一 种 关节 识别 的 方法 , 可 以 通过 

用 DTW 进行 相似 台 
台风 的 路 径 
线 , 比较 台风 
进行 台风 路 径 的 相似 性 判断 ; 
以 离 度 来 搜索 相似 路 径 的 台风 ， 这 些 文献 大 都 


[15] 在 DTW 
约束 窗 
基于 DTW 聚 
相似 水 文 时 间 序 列 
型 股市 时 间 序 列 相 
量 价 关系 问题 
一 些 运动 ， 有 效 判 
风 判 断 还 
相似 ; Ss 
生 以 及 形状 的 相似 性 


间 数 值 的 相似 ; 
文献 [21] 利 用 相 


间 序 列 , 最 具 代 表 性 的 是 动态 时 间 弯 | 


A 


(DTV)。 


三 条 三 到 4 


性 度量 方法 。 已 到了 十 征 语 音 识 另 


文献 [18] 提 


二 


要素 相似 。 本 文 研究 


等 特点 的 数 和 
为 相关 部 门 提供 准 
风 时 间 序 列 及 其 
修复 方法 ， 在 此 基础 


似 性 度量 过 程 中 各 个 要 素 重要 程度 不 同 ， 提 出 一 种 基于 加 权 动 


只 研究 台 
度 不 同 3 


断 关 键 关 节 。 目 前 ， 利 
文献 [19] 利 用 Hausdorff 距离 
化 成 平 国 


过 弯曲 时 间 轴 来 更 好 地 对 时 间 序 列 形态 进行 
别 的 一 种 主流 方 


上 县 


量 。 文献 


GhinaX 合作 期 十 
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4(U) = {S$,(D, $2(2),... S,(N)}, 
{S,(D), ,2),..., S,())} 
本 文通 过 主 成 分 分 析 法 (PCA) 计算 以 及 专家 意见 ， 确 定 


{S1(1), $2),..., Si(7)}, 
f=1,2;,07 j=12 


征 优先 匹配 原则 和 自 适应 
提出 了 一 种 
安 气 方法 ， 提 高 了 查找 
于 DTW 的 新 


度量 研究 。 


确 、 便 捷 的 辅助 决策 方法 。 


昔 述 方法 ， 然 后 给 


态 时 间 弯 
通过 实验 对 月 


E 离 的 多 元 台风 时 间 
[ 提 方 法 进行 有 效 性 分 析 。 


1 


的 多 元 台风 时 间 序 列 ， 其 长 
上 各 要 素 在 相似 性 度量 的 过 程 中 重要 程度 不 同 ， 因 此 
选择 加 权 DTW 距离 进行 相似 怕 
目标 是 设计 一 个 适合 具有 空间 性 、 


似 性 度量 方法 , 从 而 实现 


要 素 、 不 等 长 
的 相似 性 判定 ， 


先 定义 多 元 台 
风 时 间 序 的 质量 评估 与 


1 上 ， 考 虑 到 台风 时 间 序 列 的 不 等 长 性 和 相 


见 察 数据 的 时 间 点 ，( t=1,2,.…,n 
i : 观察 数据 的 变量 
c) MD : 移动 方向 
d) MV : 移动 速度 


象 的 度量 函数 


人 台风 中 心 最 低 
多 元 台风 时 间 序 列 。 一 
风 各 要 素 观测 值 8 (让 称 为 多 元 台风 时 间 序 列 ， 卓 


才 间 序列 长 度 


性 系数 阔 值 


: 一 致 性 系数 阐 值 
FP 心 近 地 面 1 min 平均 最 大 


序列 的 相似 性 度量 


个 数 ，( j=1,...,m ) 


方法 ;最 后 ， 


气压 〈 单 位 : 


风速 《〈 单 位 : 节 ) 


安 时 间 顺 序 排列 的 台 


利用 MD、MV 、P、Lo 和 La5 个 要 素描 述 台 风 ， 具 体 表示 
如 下 : 


{MD, MV,P, La,,Lo}, 
A = MD» MVsB, LasLo,), 


{MD,, MVP, La, ,Lo,} 


nr n 


定义 2 要素 依赖 。 对 于 台风 的 所 有 要 素 中 任意 两 个 要 素 
已 ， 如 果 五 已 知 ， 则 可 以 计算 出 丈 ， 则 称 到 依赖 于 五 ， 
记 作 : 态 一 五 。 台 风 等 级 G 依赖 风速 VY ， 记 作 : 了 一 G 。 


2 ”台风 时 间 序列 质量 评估 及 修复 策略 


本 章 针 对 台风 数据 中 普遍 存在 的 数据 不 完整 及 数据 不 一 致 
现象 ， 研 究 了 台风 数据 完整 性 和 一 致 性 度量 问题 ， 针 对 台风 数 
据 完整 性 及 一 致 性 修复 问题 ， 设 计 了 基于 要 素 依赖 、 近 邻 值 和 
其 他 的 台风 数据 修复 方法 。 
2.1 台风 时 间 序 列 质 量 评估 方法 
2.1.1 台风 时 间 序 列 完整 性 评估 方法 

该 台风 数据 完整 性 评估 方法 使 用 如 下 3 个 概念 : 台风 要 素 
完整 性 、 台 风 点 完整 性 、 台 风 序 列 完 整 性 。 其 定义 如 下 : 

a) 台风 要 素 完 整 性 。 台 风 要 素 完 整 性 是 指 一 个 要 素 在 时 刻 
t 时 的 完整 程度 。 对 于 A 中 任意 元 组 9 和 要 素 FF ， 要 素 值 S[F] 
的 完整 性 记 作 Crc(SLFD , 可 以 表示 为 SIF] 的 函数 ,4(S[F])， 
即 Crc(SLFD= SLF]) ， 根 据 不 同 的 应 用 ，K 可 以 具有 不 同 的 
形式 。 本 文中 函数 上 可 定义 为 
1, S[LZ] 不 为 空 值 
0, 其 他 

b) 台风 点 完整 性 。 台 风 点 完整 性 是 指 时 刻 上 时 台风 所 有 要 
素 的 完整 程度 。 对 于 A 中 任意 时 刻 上 的 数据 8 ， 8 的 完整 性 记 
作 Csc(S) 。5 的 完整 性 可 以 由 5 中 的 要 素 值 的 完整 性 来 判断 。 
则 Csc(5) 可 以 定义 为 


ztmr| 


0，5 中 主要 素 F 均 为 空 

0.2，4$ 中 主要 素 居 有 一 个 非 空 
C.(S)=10. 4，S 中 主要 素 F 有 两 个 非 空 
SC 0.6，$ 中 主要 素 F 有 三 个 非 空 
0.8，$ 中 主要 素 F 有 四 个 非 空 
1 ,中 主要 素 F 汐 非 空 


c) 台风 序列 完整 性 。 台风 序列 完整 性 是 指 一 条 台风 数据 的 
完整 程度 。 对 于 任意 一 条 台风 数据 4 ，A 的 完整 性 记 作 Cac(4) 。 
A 的 完整 性 可 以 由 A 中 时 刻 t 的 数据 5 的 完整 性 来 判断 。 则 
Cic(4) 可 以 定义 为 

Cc(A)=(NP+NPB+...+N,P)/IL 
i NN; 是 完整 度 为 B 的 数据 
5 的 个 数 。 当 Cac(4)>Q4c 时 ， 该 台风 时 间 序列 完整 性 可 修复 
的 。 
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2.1.2 台风 时 间 序 列 一 致 性 评估 方法 


该 台风 数据 一 致 性 度量 框架 使 用 如 下 2 个 概念 : 台风 点 一 


致 性 、 


由 


风 序 列 一 致 性 。 其 定义 如 下 : 


) 台风 点 一 致 性 。 人 台风 点 一 致 性 是 指 时 刻 1 时 台风 数据 中 


不 包含 语义 错误 或 相互 矛盾 的 数据 。 对 于 A 中 任意 时 刻 1 的 数 


据 $， 5 的 一 致 性 记 作 Usu(S)，Usu(5) 可 定义 为 
0，5S 中 等 级 和 风速 相互 矛盾 
1， 其 他 


Us er| 


) 台风 序列 一 致 性 。 台 风 序 列 一 致 性 是 指 一 条 台风 数据 中 


不 包含 语义 错误 或 相互 矛盾 的 数据 。A 的 一 致 性 可 以 

刻 t 的 数据 5 的 一 致 性 来 判断 。 则 Vav(4) 可 以 定义 为 
Uav(A)=N,/L 

其 中 : NN 是 满足 一 致 性 的 数量 ， 当 Var(4) > Cuv 时 ， 

间 序 列 一 致 性 可 修复 的 。 

2.2 台风 时 间 序 列 修复 方法 


A 中 


时 


风 时 


根据 文献 [22] 可 知 , 完整 性 错误 修复 结果 会 引起 一 致 性 、 时 


效 性 、 精 确 性 的 变化 ， 因 此 ， 本 文 按照 数据 完整 性 修复 、 


性 修复 的 顺序 对 台风 数据 进行 修复 。 


一 致 


台风 等 级 和 风速 两 个 要 素 之 间 存 在 依赖 关系 ， 本 考虑 要 
素 依赖 进行 完整 性 修复 ;对 于 无 法 使 用 要 素 依赖 修复 法 的 要 素 ， 


本 文 主要 考虑 风速 、 移 向 、 移 速 和 压强 四 个 要 素 ， 利 用 


之 间 的 关系 进行 完整 性 修复 ， 对 于 无 法 使 


Se 


和 风速 ， 利 用 风 压 关系 中 进行 完整 性 修复 ， 无 法 使 用 邻近 


邻近 值 修复 的 压 


值 


复 的 移 向 和 移 速 ， 根 据 经 纬度 、 观 测 值 时 间 间 隔 进行 修复 。 
台风 时 间 序 列 的 一 致 性 错误 仅 可 能 发 生 在 等 级 与 风速 之 间 ， 


本 文 利 用 要 素 依 赖 进行 一 致 性 修复 。 
具体 算法 如 下 : 


输入 : ”存在 完整 性 和 一 致 性 错误 的 台风 时 间 序 列 A(U) 


输出 : 修复 后 的 台风 时 间 序 列 
for each 4(V) do 
/根据 要 素 依赖 进行 等 级 完整 性 修复 
if( G==null&& VI=null ) 

G=G /1G 为 根据 台风 风速 确定 的 台风 等 级 


// 根 据 邻 近 值 进行 风速 、 移 向 、 移 速 、 压 强 完整 性 修复 


if( V.==null ) 


(VtV) 2, i=2,3.n, Vy < Vy 


itl 
29 


by i 
i=1 


if( MD,==null ) 


SEED | j=2,...,n, MD;, > MD ， 
MD =1| EH j=2,...,n, MD,, < MD, 
MD,, i=1 
if( MV = null ) 


MV; 1 +MV; » 
el =2,...,n, MV > MV. 


MV;_1 +MV; 
MV. 二 | -3 itl Ee 
MV,, i= 


上 


强 
用 
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if( P==null ) 


(Pi ba) 2 i=23;,n; DP <P 
i=1 


bs; 


le | ,i=2,3..n,P, >P, 


// 根 据 风 压 关 系 进行 压强 和 风速 完整 性 修复 
if( P==null,V # null ) 


P=1010_(V /16.70 
if(V = null,P # null ) 

V=6.71010-P)"™ 
// 根 据 经 纬度 进行 移 向 移 速 完整 性 修复 
if( MV = null ) 


MV = distance/ temp 
if( MD = null ) 


ALo= Lo, — Lo 
ALa= La, — La 
if(ALo>0&&aLa<0) 
angle = (90° ~ angle ) + 90° 
if(ALo<0&&aLa<0) 


angle = angle +180 
if(AaLo<0&&aLa<0) 
angle = (90 ~ angle ) +270° 
/根据 要 素 依赖 进行 一 致 性 修复 
i(V =V&&G + G,;) 
G=G// G 为 风速 WV 对 应 的 等 级 

end for 
2.3 实例 分 析 

本 文选 用 “201525” 号 台风 数据 作为 参考 台风 数据 ， 编 号 
为 1; 再 选用 “201526” 号 台风 数据 和 “201527” 号 台风 数据 ， 
编号 分 别 为 2，3，1 号 台风 原始 数据 部 分 如 下 图 所 示 : 


[{“tfbh”:”201525”, “nane”:“ 普 

琵 ", "ename” :Chanpi “is_current :1, “begin_time” :”2015-10- 

14T02:00:00”, “end_tinme” :2015-10-25T08:00:00”, “land” : [], points”: 
[{’time”:”2015-10-14T02:00:00”, “longitude” :158, 9, “latitude” :14. 0, “strong” :* 热 
带 风 暴 (T8)”, power” :8, “speed” :18, “move_dir”:“ 西 西 

北 “, “nove_speed” :22, "pressure” :998, radius?” :300, “radius10” :0, "radius12” :0, “ra 
dius?_quad’ : {ne” :300, “se” :200, “sw :240, “nw” :280}, “radius10_quad”: 

ne” :0, “se” :0, sw’ :0, nw’ :0}, “radius12_quad” | 

{"ne”:0, “se” :0, “sw” :0, nw” :0}, remark” :””, "forecast”: [{"sets”:" 中 国 ", “points’”: 


图 1 201525” 号 台风 部 分 原始 数据 


根据 2.2 节 提 出 的 方法 ， 修 复 过 程 如 下 〈 以 “201525” 号 
台风 为 例 ): 
通过 查询 ， 需 要 进行 修复 的 是 移 向 和 移 速 两 个 要 素 。 具 体 
方法 如 下 : 
a) 移 向 修复 。 台风 原始 数据 中 的 移动 方向 是 采用 十 六 风 癌 
图 记录 的 ， 为 了 便于 后 续 计算 ， 需 要 对 其 进行 数值 表示 ， 规 定 
北 为 0， 南 为 8， 从 北 顺 时 针 到 南 依次 是 1 一 7; 从 北 逆 时 针 到 
南 依次 是 9 一 15。 
根据 经 纬度 计算 出 角度 ， 对 照 十 六 风向 图 以 及 数值 转换 关 
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系 补充 移 向 。 例如 “2015-10-14T14:00:00” 时 刻 , 经 纬度 分 别 为 
156.9、15.0， 后 一 时 刻 的 经 纬度 分 别 为 154.4、15.8， 根 据 2.2 
节 方法 计算 出 角度 为 288.3294336171956, 对 照 十 六 风向 图 得 出 
移 向 是 西西 北 ， 数 值 化 结果 为 11。 

b) 移 速 修复 。 根 据 经 纬度 以 及 时 间 间 隔 计算 移 速 。 例 如 
“2015-10-14T14:00:00” 时 刻 ， 经 纬度 分 别 为 156.9、15.0， 后 
一 时 刻 “2015-10-15T02:00:00” 其 经 纬度 分 别 为 154.4、15.8， 
计算 出 移 速 取 整 为 28。 

利用 本 文 提 出 的 修复 方法 处 理 之 后 的 数据 如 下 : 


1 号: {{11.0,22.0.998.0,14.0,158.9}1, 
让 i 
{11.0,23.0,990.0,15.5,152.6}s,..., 
名 .0.00982 .031.1.158.0je 
号 : {{11.0,20.0,998.0,4.6,159.7}, 
{11.0,20.0,998.0,5.1.158.8}»...., 
{11.0,20.0,985.0,7.7,153.4}s,..., 
人 .0.40.0.998.0.21 3} 
号 : {{11.0,25.0.998.0,9.1,138.0}1, 
{11.0; oo os 723. ee 
{11.0,19.0,990.0,11.2,132.3 }s,. 
ee 1g 8 


图 2 修复 后 台风 数据 


其 中 , 每 个 元 素 的 下 标 表示 该 点 在 整个 时 间 序 列 中 的 位 置 ， 
最 后 一 个 下 标 值 表示 整个 时 间 序 列 的 长 度 。 

经 检验 , 移 向 的 均 方 根 误差 (root mean square error, RMSE ) 
值 为 0.187867287325544, 移 速 的 RMSE 值 为 0.264797369811 ， 
该 修复 方法 精确 性 高 。 


3 ”基于 加 权 DTW 距离 的 台风 时 间 序 列 的 相似 性 度量 


台风 具有 空间 性 、 季 节 性 和 多 要 素 等 特点 。 本 文 综合 考虑 
移 向 、 移 速 、 压 强 、 经 度 和 纬度 五 个 要 素 ， 针 对 各 要 素 重 要 程 
度 不 同 且 台风 时 间 序 列 不 等 长 的 特点 ， 本 文采 用 加 权 DTW 距 
离 进行 相似 性 度量 。 

WO 
素 权 值 设计 ， 加 权 DTW 距离 计算 。 季 节 相 似 容易 实现 ， 不 再 
效 述 ， 本 章节 主要 介绍 台 eh 
算 两 个 方面 。 
3.1 台风 要 素 权重 设计 

本 文 利用 主 成 分 分 析 法 ， 确 定 描 述 台 风 的 五 个 要 素 ， 分 别 
是 台风 的 移动 方向 、 移 动 速度 、 压 强 、 纬 度 和 经 度 ， 它 们 的 权 
重 分 别 用 WW，WW，W，W 和 WW 表示。 权重 计算 方法 是 利用 
PCA 方法 计算 出 相关 系数 矩阵 及 其 特征 向 量 和 特征 值 , 根据 设 
定 的 主 成 分 贡献 率 阔 值 ， 判 断 主 成 分 的 个 数 ， 再 根据 主 成 分 中 
各 要 素 的 系数 计算 出 台风 各 个 要 素 的 权重 。 计 算 过 程 如 下 : 

a) 对 原始 台风 数据 进行 标准 化 处 理 

b) 计算 相关 系数 矩阵 

c) 计算 相关 系数 矩阵 的 特征 向 量 和 特征 值 

特征 向 量 B 可 以 表示 成 如 下 所 示 : 


du CD 1 du is 
a dd» Cs (02 dss 
B=| qs dy ds dx dss 
da ds da Gu dss 
Cs ds ds ds dss 


,0 人。 


黄 冬 梅 ， 等 : 一 种 多 lo 

将 特征 值 按照 从 大 到 小 的 顺序 排列 ， 依 次 是 和, 帮 …,4， 
d) 选择 主 成 分 并 计算 其 贡献 率 

主 成 分 表示 如 下 所 示 : 


M, =aiMP+aMV +...+ ALo 
M, =aMDP+a MVv +..+ Lo 


M, =a,MD+a,MV +..+a,Lo 


每 个 主 成 分 中 台风 要 素 的 系数 是 相关 系数 矩阵 特征 向 量 的 
主 成 分 的 贡献 率 如 下 所 示 : 


C=4/(N + +...+h) 
Ca = 邦 /(4 + 石上 + 十 力 ) 


C =4/0+1++4) 
e) 计算 台风 各 个 要 素 的 权 值 
假设 主 成 分 有 n(n<5) 个 ， 则 移 向 、 移 速 、 压 强 、 纬 
和 经 度 的 权 值 分 别 如 下 : 


CI 
人 了 |a， | GG Cc 


列 


洱 


平 la, © 
Ga C++ C+C,+...+C, 
W,= a 有 G jp 
~ altlaslt:+tlas) C+ t+tO, 
9 la x CG 
lt ae C+C,+...+C, 
W, = la CI 
laul+tlaslt.+la, G6 F.C 
Cm C: 


十 X 
lasltlaslt..+la,, 和 全 二 全 Go 


3.2 ”加权 DTW 距离 计算 
对 于 不 同 的 台风 ， 人 台风 各 要 素 对 相似 性 度量 的 影响 不 同 ， 
本 文 利用 加 权 DTW 距离 进行 相似 性 度量 ， 具 体 定义 如 下 ; 
设 时 间 序 列 X={S,(D, 5,(2, 5,G), 5,(4), 5,(5))} ， 


Y={5,(D，5,(2)，S,(3)，S,(D，5,(5)} ， 则 入 ， 了 的 加 权 DTW 


距离 定义 见 式 (1)。 

0 i=0,j=0 

oo i=0 j#0orj=0, i*#0 (1) 
d (x,y))+ Min else 


Min=min{7 (i—1,j—1),r(i-1,)7),7(i,7 -1)} 


r(i,))= 


其 中 :对 于 i=1,2,...,m j=42,..sn ， 用 r(i,j) 表示 
r(X(1:),7(1: 7)) ，4d(%,y,) 表 未 和 yy) 之 间 的 基 距 离 ， 可 以 
根据 情况 选择 不 同 的 距离 度量 ， 本 文 利 用 如 下 公式 进行 计算 ， 


et 


+W (BP, +W (La 1a) +W (Lo -10)) 


多 元 台风 时 间 序 列 相似 性 度量 的 算法 如 下 : 
输入 : 历史 台风 数据 arr2 和 参考 台风 数据 arrl; 
输出 :历史 台风 与 参考 台风 的 距离 d; 
def Distance(i,j,arr1,arr2, W): 
if(i==0 and j ==0): 

return 0 


a 
录用 稿 黄 冬 梅 ， 等 : 一 种 多 Ghina yi, 
elseif (i==0 or j==0): 括 台 风 的 中 英文 名 称 、 发 生 时 间 、 结 束 时 间 、 经 纬度 、 强 度 、 
return float("inf"') 风速 、 移 动 方 向 、 移 动 速度 、 压 强 等 要 素 的 值 。 
else: 4.2 实验 结果 与 分 析 
d=0 相似 台风 是 指 两 条 台风 季节 相似 及 满足 加 权 DTW 距离 要 
datal =arrl[j] 求 。 随 机 选择 一 条 台风 数据 为 参考 数据 ， 编 号 为 1， 再 将 收集 
data2 =arr2[j] 到 的 历史 台风 数据 分 别 编号 为 2，3.……. ，N (WN 表示 历史 台 
fori iin range(len(arrl[0])): 风 数 据 的 数量 )。 如 果 编 号 2,3...N 都 满足 季节 相似 ， 则 计算 其 
d=d+pow(datal[i]-data2[i],2)*WI[i] 与 参考 数据 之 间 的 距离 ， 记 为 dz2，4i3..….……. ，diw ， 如 果 距 离 
d=powf(result,0.5) 越 小 ， 则 两 条 台风 越 相似 。 本 文 设 计 了 3 个 实验 ， 实 验 一 台风 
return result+ 时 间 序 列 质 量 评估 与 修复 ， 实 验 二 是 多 元 台风 时 间 序 列 相似 性 
min(Distance(i-1,j-1,arrl,arr2,W) 度量 方法 的 有 效 性 验证 ; 实验 三 是 多 元 台风 时 间 序 列 相似 性 度 
Distance(i-l,j,arrl,arr2,W)， 量 方法 的 实用 性 证 明 。 
Distance(i,j-1,arrl,arr2, W)) 实验 1 台风 时 间 序 列 质 量 评估 与 修复 
算法 分 析 : 设 参考 台风 数据 长 度 为 WMW， 历史 台风 数据 长 度 本 实验 包含 两 组 实验 ， 第 一 组 实验 是 对 需要 进行 实验 的 数 
为 N， 该 算法 的 时 间 复 杂 度 为 O(MN) 。 据 进行 完整 性 评估 及 修复 ， 第 二 组 实验 是 对 本 文 提 出 的 台风 时 
3.3 ”实例 分 析 间 序 列 修复 方法 进行 有 效 性 验证 。 
本 节 使 用 的 数据 是 2.3 节 中 处 理 之 后 的 数据 ， 具 体 计算 过 1) 台风 时 间 序 列 完 整 性 评估 与 修复 
程 如 下 所 示 (以 编号 1 和 2 的 数据 为 例 ): (1) 完整 性 闵 值 确定 
首先 计算 1 号 台风 各 个 要 素 的 权 值 ， 本 文采 用 主 成 分 分 析 实验 数据 采用 “195526”“200008”、“200917”、“201117” 
法 ， 得 到 的 结果 如 表 1 所 示 。 四 条 台风 数据 。 利 用 本 文 提出 的 方法 进行 完整 性 评估 ，“200008” 
表 1 主 成 分 、 贡 献 率 及 权 值 结 号 (编号 3) 台风 的 完整 性 是 0.6,“200917” (编号 4) 号 台 
M, M, M, M, M; WwW 风 的 完整 性 约 为 0.8， 将 “195526” 号 〈 编 号 2) 台风 数据 进行 
MD 0.1291 -0.1905 0.6530 0.1614 0.7033 0.0803 区 改 ， 使 其 完整 性 为 0.4, 分 别 计算 三 条 台风 与 “201117” 号 ( 编 
MV 0.4315 0.8247 0.3286 -0.0677 -0.1454 0.2876 号 1) 台风 之 间 的 DTW 距离 ， 如 表 2 所 示 。 
P 0.8544 -0.3985 -0.1897 -0.2729 -0.0260 0.4426 表 2 考虑 路 径 相 似 的 距离 (“201117”) 
La -0.0526 0.3369 -0.5441 -0.3441 0.6851 0.0594 编号 加 权 DTW 距离 
Lo 0.2535 0.1064 -0.3655 0.8812 0.1194 0.1301 1,2 37806.25047811353 
C 0.7737 0.1824 0.027 0.0103 0.0066 1.3 45169.845710044625 


1,4 77485.42226788626 


然后 计算 1 号 台风 和 2 号 台风 之 间 的 距离 为 


表 2 可 知 ， 与 “201117” 号 台风 相似 度 从 大 到 小 的 排列 
i 2 ，\2 是 “195526” 号 、“200008” 号 和 “200917” 号 , 而 已 知 “200008” 

+W (Ps -Ps) + 有 (La -La ss) +W (Lom -Los) 2、 i A 症 、 口 粮 
0.0803x (2.0—2.0) +0.2876x(42.0-40.0 + 号 台风 与 “201117” 号 台风 更 相似 ， 所 以 存在 错误 ， 即 在 数据 
补充 过 程 中 ,完整 性 低 于 0.6, 修复 的 数据 错误 率 较 高 。 故 本 文 
中 ， 只 有 当 完 整 性 大 于 等 于 0.6 时 ， 修 复数 据 进 行 相似 性 度量 


W (MD, -MD W (MV, -MTV . 
= + 全 
d (xe,y5)= ( 到 5) 人 67 .| 


= |0.4426x(982.0—998.0) +0.0594x(31.1-21.7) 
+0.0301x(158.0—137.7) 


将 该 结果 代入 式 (2), 得 出 结果 为 578.4718754376307, 同 ” 才 更 准确 。 
理 ， 计 算出 1 号 台风 和 3 号 台风 的 距离 为 707.0062487978438， (2) 台风 时 间 序 列 完整 性 评估 
于 1 号 台风 和 2 号 台风 之 间 的 距离 更 小 ， 所 以 这 两 条 台风 更 本 实验 选取 的 数据 是 下 面 两 个 实验 需要 使 用 的 数据 ， 分 别 
相似 。 是 “201117” 号 台风 、“200008” 号 台风 、“195526” 号 台风 、 


“199327 ”号 台风 “200917? 号 台风 “201323 ”号 台风 “197010?” 


号 台风 、“201312” 号 台风 、“200713” 号 台风 、“200414” 号 台 


4 ”实验 结果 与 分 析 


4.1 实验 环境 与 实验 数据 风 、“195615” 号 台风 、“199112” 号 台风 、“195316” 号 台风 、 

实验 环境 为 MATLAB R2013a，IntelliJ IDEA 2016.3.4.lnk， “198510” 号 台风 和 “199111” 号 台风 15 条 台风 数据 。 这 些 数 
Win7SP1，1TB 硬盘 ，8 GB 安装 内 存 ，IntelG) CoreIM i7-3770 ”” 据 主要 存在 完整 性 问题 ， 修 复 前 后 台风 数据 完整 性 对 比 图 如 图 
CPU。 3 所 示 。 


选取 1945 一 2016 年 的 台风 数据 作为 研究 对 象 。 台 风 数 据 包 
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图 3 修复 前 后 台风 时 间 序列 完整 性 对 比 图 5 风 压 修复 方法 精确 性 


上 图 可 知 , 修复 之 后 台风 时 间 序 列 的 完整 性 是 100%, 便 
于 后 续 台 风 相 似 性 度量 ， 提 高 相似 性 度量 的 准确 率 。 

2) 台 风 时 间 序列 修复 方法 精确 性 验证 

实验 数据 采用 “200713”“200817”、“200917”、“201011”、 
“201117”“201215”“201312”“201323”“201416” 和 “201521?” 
10 条 台风 数据 。 

(1) 完整 性 修复 方法 精确 性 验证 

根据 文献 [26] 可 知 ， 修 复 精 确 性 可 用 如 下 公式 度量 


其 中 : e 是 台风 要 素 原 有 的 真实 值 ，E& 是 修补 的 值 ，m 是 台 反 
的 长 度 ，RMSE 的 值 越 小 ， 说 明 精确 性 越 高 ， 反 之 ， 精 确 性 越 


图 6 经 纬 计算 修复 方法 精确 性 


低 。 (3) 一 致 性 修复 方法 精确 性 验证 
将 原来 完整 的 数据 随机 抽取 10%， 用 null 值 代替， 利用 本 随机 抽取 原来 正确 数据 的 10%， 将 等 级 数据 改 成 与 风速 不 


文 提 出 的 方法 进行 修复 ， 采 用 邻近 值 修复 修复 方法 的 精确 性 如 一致 的 数据 ， 得 到 的 一 致 性 修复 方法 的 精确 性 如 图 4 所 示 。 
图 4 所 示 。 


250 -5 
0.8 
200 
0.7 
0.6 家 fi 
05 9 a ld i 
= 等 级 曙 ioo ， 上 修复 前 一 至 数量 
2 0, 卓 速 度 目 修 复 后 一 致 数量 
别 目 移 向 0 了 目 总 长 度 
虽 移 如 
0.1 j 个 - 寺 
I 和 县 呈 基 中 目 压 强 


AAA > 
y ee Ee gr > Hyg > y RY sp 
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台风 编号 
图 4 邻近 值 修复 方法 精确 性 图 7 一 致 性 修复 方法 的 精确 性 


将 原来 完整 的 数据 随机 抽取 10%， 抽 取 的 数据 每 组 至 少 有 图 4 可 知 ， 每 个 要 素 的 RMSE 的 值 都 低 于 0.8; 图 5 可 
两 个 null 值 是 相 邻 的 ， 将 其 中 的 压强 和 风速 用 null 值 代替 ， 利 。 知 ， 风 速 和 压强 的 RMSE 的 值 均 小 于 4， 图 6 可 知 ， 移 向 移 速 


用 风 - 压 关系 修复 方法 进行 修复 ， 其 修复 精确 性 如 图 5 所 示 。 的 RMSE 的 值 均 小 于 0.7， 所 以 ， 本 文 提出 的 完整 性 修复 方法 
将 原来 完整 的 数据 随机 抽取 10%， 抽 取 的 数据 每 组 至 少 有 修复 方法 精确 性 较 高 。 

两 个 null 值 是 相 邻 的 ， 将 其 中 的 移 向 移 速 用 null 值 代 奉 ， 利 用 实验 2 多 元 台风 时 间 序 列 相似 性 度量 方法 有 效 性 验证 

经 纬 计算 修复 方法 进行 修复 ， 其 修复 精确 性 如 图 6 所 示 。 本 实验 包含 两 组 实验 ， 第 一 组 数据 是 “201117” 号 台风 数 


据 及 其 相似 路 径 台风 数据 ， 第 二 组 数据 是 “201323 ”号 台风 数 
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据 及 其 相似 路 径 人 台风 数据 。 利 用 本 文 的 方法 ,只 考虑 路 径 相似 ， 
即 移 向 移 速 两 个 要 素 ， 计 算 加 权 DTW 距离 。 

1) 与 “201117” 号 台风 路 径 相 似 的 台风 

本 实验 设置 贡献 率 的 阔 值 为 90%， 权 值 计算 结果 如 表 3 前 
两 列 所 示 。 


表 3 201117” 号 与 “201323” 号 主 成 分 、 贡 献 率 及 权 值 结果 


主 成 分 ” 权 值 ” 主 成 分 ” 权 值 
移 向 (MD) 0.0100 0.0099 0.2257 0.1881 
移 速 (MV) 0.9999 0.9901 0.9742 0.8119 
贡献 率 CC ) ”0.9967 0.9212 


选择 “201117” 号 台风 数据 为 参考 数据 , 编号 为 1， 由 文献 
[26] 可 知 , 与 “201117” 号 台风 相似 的 台风 按照 相似 度 从 高 到 低 
依次 是 “200008” 号 , “195526” 号 ,“199327” 号 和 “200917” 
号 ， 分 别 编号 为 2，3，4，5。 只 考虑 路 径 相 似 ， 由 表 3 可 知 ， 
移 向 移 速 的 权 值 分 别 为 0.0099 和 0.9901, 分 别 计算 出 台风 2 与 
台风 1 的 距离 ， 台 风 3 与 台风 1 的 距离 ， 台 风 4 与 台风 1 的 距 
离 ， 台 风 5 与 台风 1 的 距离 ， 计 算 结 果 如 表 4 所 示 。 


表 4 考虑 路 径 相似 的 距离 (“201117”) 
编号 加 权 DTW 距离 
1,2 206.98736773989432 


1,3 211.7052704239524 
1,4 317.9152515180596 
1,5 849.5429360222836 


根据 计算 结果 可 知 , 距离 大 小 按照 从 低 到 高 排序 为 1,2、1,3、 
1,4、1,5， 即 与 “201117” 号 台风 相似 台风 按照 相似 性 从 高 到 低 
分 别 为 “200008” 号 , “195526” 号 , “199327” 号 和 “200917” 
号 ， 与 文献 [26] 给 出 的 结果 相同 。 

2) 与 “201323” 号 台风 相似 的 台风 

选择 “201323” 号 台风 数据 为 参考 数据 , 编号 为 1， 根据 文 
献 [20], 与 “201323” 号 台风 相似 的 台风 按照 相似 度 从 高 到 低 依 
次 是 “197010” 号 ,“201312” 号 和 “200713” 号 ， 分 别 编号 为 
2，3，4， 只 考虑 路 径 相 似 性 ， 由 表 3 后 两 列 可 知 ， 移 向 移 速 的 
权 值 分 别 为 0.1881 和 0.8119， 加 权 DTW 距离 计算 结果 如 表 5 
所 示 。 


表 5 考虑 路 径 相 似 的 距离 
编号 加 权 DTW 距离 
1.2 328.82862494660293 


1,3 445.0503741631772 
1,4 613.0328465781108 


根据 计算 结果 可 知 , 距离 大 小 按照 从 低 到 高 排序 为 1,2、1,3、 
14， 即 与 “201323” 号 台风 相似 台风 按照 相似 性 从 高 到 低 分 别 
为 “197010” 号 ,“201312” 号 和 “200713” 号 ， 与 文献 [20] 给 
出 的 结果 相同 。 
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上 述 两 个 实验 表明 ， 本 方法 能 够 有 效 判断 台风 路 径 相 似 。 
以 上 两 篇 文章 中 的 方法 均 没 有 考虑 台风 的 季节 性 特点 ， 根 据 本 
文 的 方法 , 考虑 季节 相似 , 与 “201117” 号 相似 的 只 有 “195526” 
号 和 “200917” 号 , 与 “201323” 号 台风 相似 的 只 有 “200713” 


‘a 


实验 3 多 元 台风 时 间 序 列 相似 性 度量 方法 实用 性 验证 

选择 “200414” 号 台风 数据 为 参考 数据 , 编号 为 1， 本 实验 
设置 贡献 率 的 六 值 为 90%， 利 用 2.1 节 权 重 计算 方法 ， 计 算出 
移 向 移 速 的 权 值 分 别 为 0.0315, 0.9685。 如果 只 考虑 路 径 相似 ， 
得 到 的 5 个 相似 台风 按照 相似 度 由 大 到 小 排序 分 别 是 “195615” 
号 、“199112” 号 、“195316” 号 、“198510” 号 和 “199111” 号 ， 
将 其 分 别 编 号 为 2、3、4、5、6， 加 权 DTW 距离 如 表 6 所 示 。 


表 6 考虑 路 径 相似 的 距离 (“200414””) 
编号 加 权 DTW 距离 
1,2 179.7088665502276 


1,3 195.154657468981 

1,4 196.21462224736027 
1,5 200.0512292651246 
1,6 201.13529612422815 


本 文 除了 考虑 路 径 ， 考 虑 到 台风 的 空间 性 以 及 台风 灾害 的 
影响 ， 还 综合 考虑 了 压强 、 经 纬度 三 个 要 素 ， 能 够 更 加 全 面 的 
度量 相似 台风 ， 进 而 有 效 预 判 台风 带 来 的 灾害 程度 ， 提 前 在 相 
关 地 区 做 好 防护 措施 。 移 向 、 移 速 、 压 强 、 纬 度 和 经 度 的 权 值 
分 别 为 0.0008、0.0052、0.9926、0.0012 和 0.0002， 计 算 的 加 权 
DTW 距离 结果 如 表 7 所 示 。 

表 7 综合 多 要 素 的 距离 
编号 加 权 DTW 距离 
1,2 1880.052322218788 


1,3 1259.9505176627304 
1,4 1840.153387085322 
1;3 1099.7791635028807 
1,6 1183.6785528177586 


上 表 可 知 ， 与 “200414” 号 台风 相似 的 5 条 台风 相似 度 
大 到 小 排序 分 别 是 “198510” 号 、“199111” 号 、“199112” 
号 “195316” 号 和 “195615” 号 ， 而 与 “200414” 号 台风 路 径 
最 相似 的 是 “195615” 号 台风 。 由 下 图 中 国 台风 网 中 各 个 台风 
的 路 径 及 强度 信息 可 以 看 出 ,，“195615” 号 台风 强度 变化 与 
“200414” 号 台风 最 不 相似 ,“198510” 号 台风 强度 变化 与 
“200414” 号 最 相似 。 

根据 文献 [27] 可 知 ,“200414” 号 台风 云 娜 在 浙江 造成 的 灾 
害 主 要 表现 在 以 下 三 个 方面 : a) 大 风 导 致 大 量 房屋 倒塌 ， 造 成 
的 死亡 人 数 占 总 死亡 人 数 的 66. 5%; 人 员 被 刮 倒 而 死亡 的 人 数 
5 总 死亡 人 数 的 5. 5%; 部 分 电 杆 吹 倒 压 死 或 因 触 电 而 死 的 人 
数 占 总 死亡 人 数 的 3%;b) 由 于 台风 的 强风 和 低 气 压 的 作用 , 使 
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海水 向 海岸 方向 强力 堆积 ， 潮 位 猛 涨 ， 造 成 沿海 地 区 堤防 损坏 ” 4.3.2 参数 性 能 对 比 
4059 处 563 km, 堤防 决口 1222 处 88 km;c) 台风 强 降水 使 44. 本 文 对 几 种 相似 台风 的 判断 方法 进行 考虑 要 素 、 正 确 性 以 
4 万 名 群众 被 洪水 围困 ， 还 诱发 了 强 地 质 灾害 。 文 献 [28] 指 出 ， ”及 效率 方面 的 比较 ， 得 到 的 结果 如 下 表 所 示 ， 由 表 中 结果 可 以 
“198510” 号 台风 造成 了 广西 持续 一 周 的 暴雨 一 一 大 暴雨 的 天 看 出 ， 本 文 的 方法 在 各 个 方面 都 优 于 其 他 几 种 方法 。 

气 ; 文献 [29] 指 出 ,，“198510” 号 台风 引发 了 泥石流 ， 破 坏 了 一 表 8 参数 性 能 比较 
批 水 利 工程 ,可 见 ,“198510” 号 台风 带 来 的 灾害 程度 与 “200414” 雍 称 相似 。 沽 雇 灾 。 党 囊 兴 ” 沽 下 地 
号 合 风 相 似 。 性 度 最。 害 各 度 。 地 影响 再 位置 开 
基于 WebGIS 的 分 
是 否 否 是 较 正确 高 
析 系 统 [19] 
相似 离 度 原理 [21] 是 否 否 否 正确 较 高 
基于 ArcGIS 的 筛 
是 是 否 否 正确 高 
选 方法 [24] 
本 文 方法 是 是 是 是 正确 高 
图 8 台风 信息 5 ”结束 语 

文献 [1] 表 明 ， 超 强 台 风 “ 英 兰 蒂 ” 由 于 其 超 强 强度 ， 致 使 本 文 首先 结合 台风 的 特点 ， 利 用 多 元 时 间 序 列 对 台风 数据 
福建 、 浙江、 江西、 上 海 、 江 苏 5 个 省 (市)375.5 万 人 受灾 ， 进行 了 形式 化 定义 ， 实 现 了 具有 空间 性 、 季 节 性 以 及 多 要 素 特 
44 人 死亡 失踪 ， 直 接 经 济 损 失 316.5 亿 元 。 由 此 可 见 ， 对 于 台 ”点 的 数据 的 定量 描述 ， 对 台风 数据 中 存在 的 完整 性 和 精确 性 问 
风灾 害 防护 方面 的 应 用 ， 不 仅 需 要 考虑 台风 路 径 ， 还 需要 考虑 。 题 进行 评估 和 修复 ， 然 后 以 其 为 基础 ， 提 出 了 一 种 基于 加 权 
其 强度 ， 因 此 ， 直 接 比 较 路 径 相似 性 是 不 全 面 的 ， 需 要 综合 考 ” DTW 上 距离 的 多 元 台风 时 间 序 列 的 相似 性 度量 方法 ， 实 验 结果 
虚 多 个 因素 ， 实 现 台风 灾害 的 准确 预测 ， 提 前 做 好 防护 措施 。 表明 ， 本 文 提出 的 台风 数据 质量 评估 方法 能 够 提高 相似 性 度量 
4.3 性 能 分 析 的 正确 率 和 灾害 预测 的 正确 性 ， 相 似 性 度量 方法 能 够 根据 不 同 
4.3.1 时 间 性 能 对 比 的 参考 台风 设置 不 同 的 权 值 以 适合 各 种 相似 台风 的 判断 ， 能 够 

选择 “194930” 号 、“195820” 号 、“196911” 号 、“197613” 有 效 判断 台风 造成 的 灾害 。 下 一 步 工作 是 根据 历史 台风 数据 进 
号 “198510”“199215” 号 “200414” 号 “201117” 号 “201323” 行 实时 台风 走势 、 强 度 等 预测 ， 为 相关 部 门 提供 防 灾 减 灾 、 辅 
号 和 “201521” 号 十 条 数据 进行 运行 时 间 的 对 比 ， 利 用 本 文 的 。” 助 决策 。 

a ,分 别 考虑 移 向 \ 移 速 2 个 要 素 和 在 此 基础 上 ， 参考 文献 

增加 考虑 压强 、 经 纬度 共 5 个 要 素 所 需要 的 相似 性 度量 时 间 ， 

得 到 的 结果 如 图 9 所 示 。 1] 中 国 气象 局 .2016 年 中 国 气 候 公报 . [EB/OL]. [2017-09-04]. http://m. cma. 
gov. cn/root7/auto13139/201705/t20170525_415102. html. 

. i 2] 李 建 中 , 王 宏 志 , 高宏, 大 数据 可 用 性 的 研究 进展 []]， 软件 学 报 ,2016， 

6 | “考虑 5 个 要 素 27 (7): 1605-1625. 
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图 9 ”时间 性 能 对 比 图 
图 9 可 以 看 出 ,增加 考虑 压强 、 经 度 和 纬度 三 个 要 素 后 ， 
运行 时 间 增 加 ， ee 实验 表明 , 本 文 的 方法 ， 
考虑 要 素 更 加 全 面 ， 且 运行 效率 几乎 不 变 。 
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